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Abstract 

We study the physics capabilities of the NuMI beamline with an off-axis highly- 
segmented iron scintillator detector and with the inclusion of the currently under study 
proton driver upgrade. We focus on the prospects for the experimental determination 
of the remaining neutrino oscillation parameters, assuming different outcomes for ex- 
periments under way or in preparation. An optimization of the beam conditions and 
detector location for the detection of the i/„ — * ^-transitions is discussed. Different 
physics scenarios were considered, depending on the actual solution of the solar neu- 
trino puzzle. If KamLAND measures Am^, we find it possible to measure both |C/ e 3| 2 
and the CP violating phase 5 within a viable exposure time, assuming a realistic detec- 
tor and a complete data analysis. Exposure to both neutrino and antineutrino beams 
is necessary. We can, in addition, shed light on Am|> if its value is at the upper limit 
of KamLAND sensitivity (i.e., the precise value of Amq remains unknown even after 
KamLAND). If the solar neutrino solution is not in the LMA region, we can measure 
|[7e3| 2 and determine the neutrino mass hierarchy. The existence of the proton driver 
is vital for the feasibility of most of these measurements. 



1 Introduction 



There is hope that the reconstruction of the neutrino mass matrix will shed light on some of 
the most relevant open questions faced by high energy physics today, including the origin of 
the neutrino mass, the flavor-puzzle, and the origin of the asymmetry between matter and 
antimatter, among others. 

In order to start tackling these issues, "precision measurements" of the neutrino parame- 
ters are required: simply knowing that the neutrino masses are tiny is not enough. Knowing 
that there are three non degenerate massive neutrinos is also not enough to prove that there 
is CP violation in the lepton sector, although it would be rather surprising if it turned out 
otherwise. Forthcoming neutrino experiments, therefore, have to be designed not simply to 
test the oscillation hypothesis, but also to perform precise measurements of the oscillation 
parameters. 

We currently know that there are (at least) three neutrino "flavors" - the electron, 
the muon and the tau neutrinos - and the current neutrino oscillation data strongly sug- 
gest that the flavor eigenstates differ from the mass eigenstates and therefore the neutrinos 
can mix. Furthermore, the atmospheric, solar and reactor data point to two hierarchically 
different mass-squared differences plus a small "connecting angle," such that, to good ap- 
proximation, both the atmospheric and solar neutrino puzzles can be solved by assuming 

<r-> z/ T -oscillations and v e z/ ot her-oscillations, respectively. Forthcoming results from so- 
lar, atmospheric, accelerator, and reactor neutrino oscillation experiments will measure the 
size of the mass-squared difference and mixing angle that drive the solar neutrino oscillations 
with good precision and the mass-squared difference and the mixing angle involved in the 
atmospheric oscillation to about 10% [l|] (and there should be direct confirmation that muon 
neutrinos oscillate to tau neutrinos by direct observation of a tau appearance 0). Further- 
more, we will also definitively open (or close) the door to extra new physics in the neutrino 
sector by confirming (or not) the LSND anomaly || with MiniBooNE data [[|]|] 

In addition, some non-oscillation neutrino experiments are aiming at determining the 
absolute value of the neutrino masses using direct searches || and the nature of the neutrino 
mass: are they Majorana or Dirac particles? ||. The knowledge of the neutrino mass is 
important to decide whether neutrinos contribute significantly to the total energy of the 
Universe. 

In view of the above results and the ongoing experimental programs, the ultimate goal for 
the next generation of neutrino experiments should be to test CP violation in the neutrino 
sector, if the parameters of the solar neutrino solution allow such a determination. One way 
to achieve this goal is to devise experiments that: 

• are sensitive to the sub-dominant — > v e channel for the atmospheric L/E v , and can 
measure \U e ^\ from both — > v e and — > v e transitions; 

• are capable of determining the neutrino mass pattern from matter effects; 
*We will not consider the LSND anomaly in this paper. 
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• ultimately, can test CP invariance in the leptonic sector. 

Such experiments require, in general, relatively intense and mono-energetic beams. As 
first suggested in J7|, this could be achieved with off-axis neutrino beams. 

The literature contains a significant amount of work || || which discusses the potential 
of generic "super-beams," with neutrino energies ranging from 1 to 50 GeV and baselines 
spanning from 200 to 7000 kilometers. Low energy neutrino beams have also been studied 
TI]| ]. In this work, we explore which of the questions above can be addressed by using the 



NuMI beam line with a new proton driver and the construction of an off-axis detector. We 
will concentrate exclusively on what one may hope to learn by looking for — *> u e appear- 
ance in a highly segmented iron-scintillator detector. We will conclude, in agreement with 
some of these previous works, that intense beams can significantly improve our knowledge 
of the neutrino oscillation parameters, including (depending on the solar solution) some sen- 
sitivity to a CP violating phase. However, the ultimate sensitivity to some of the neutrino 
parameters, in particular the CP violating phase, will require the purity and intensity of 
neutrino factory beams |]TT] . 



The paper is organized as follows. First, we briefly review the neutrino mixing matrix 
and oscillation probabilities, and describe the current knowledge of neutrino masses and 
mixing angles. In Sec. 3, we describe the off-axis neutrino beam, and discuss where an off- 
axis detector should be located in order to maximize its physics capabilities. In Sec. 4, we 
describe in detail the detector we will be considering, and discuss reconstruction efficiencies 
and strategies for reducing the number of background events. In Sec. 5, we discuss the physics 
capabilities of such a setup, for different values of the solar mass-squared difference. We find 
that while some information can be obtained from a "neutrino" beam, it is imperative to 
run with an "antineutrino" beam as well. Furthermore, in order to extract more "exciting" 
physics (such as CP-violation) out of the experiment, a new intense proton source and a 
large detector are required. In Sec. 6 we summarize our results and compare our findings 
with similar studies at different beamlines. 



2 Neutrino Mixing and Oscillations 

The presence of non-zero masses for the light neutrinos introduces a leptonic mixing matrix, 
U, analogous to the well known CKM quark mixing matrix, and which in general is not 
expected to be diagonal. The matrix U connects the neutrino flavor eigenstates with the 
mass eigenstates: 

W) = ££yi>i>, (2.i) 

i 

where a denotes the active neutrino flavors, e, /i or r, while % runs over the mass eigenstates. 
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It is "traditional" to define the mixing angles #12,13,23 in the following way: 
i ;u! i/ ]2 = . tan 2 #23 

Well 

while 
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^(U^U^U^U*^ = sin #12 cos #12 sin #23 cos #23 sin 9 13 cos 2 6> 13 sin 5, (2.3) 

defines the CP-odd phase 5. For Majorana neutrinos, U contains two further multiplicative 
phase factors, but these are invisible to oscillation phenomena. 

In order to relate the mixing angles and mass-squared differences to the parameters 
constrained by experiments, it is convenient to define the neutrino masses such that m 2 < m 2 
with Am 2 2 < I Am 2 3 23 |[] (the data, in fact, point to Am 2 2 -C | Am 2 3 23 |). With this definition, 



the "solar angle" 9 Q ~ #12, while the "atmospheric angle" ^atm 
experiments constrain |£/ e3 | 2 . The solar mass-squared difference Am^ 



#23- Furthermore, reactor 
Am 2 2 , while the 



|Am 2 3 | 



|Am 2 3 | 



It is important to note 



atmospheric mass-squared difference is Am 2 tm 
that m 3 can be either larger or smaller than m 2 , mf. 

The oscillation probability P{y a — > vp) is given by the absolute square of the overlap of 
the observed flavor state, \vp)-, with the time-evolved initially-produced flavor state, \v a ). In 
vacuum, it yields the well-known result: 
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(2.6) 



P(f Q -> Z//3) = P(z//3 -» f a ) = P C P-cven(^ ~> f/j) ~ PcP-odd(^a ~> ^) 

where, by CPT invariance, P[y a — > vp) = P{pp — > v a ). In vacuum the CP-even and CP-odd 
contributions are even and odd, respectively, under time reversal: a *->■ [3. 

If the neutrinos propagate in ordinary matter, these expressions are modified. The prop- 
agation of neutrinos through matter is very well described by the evolution equation 

,2 



dt 
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(2.7) 
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where A/(2E U ) is the amplitude for coherent forward charged- current scattering of v e on 
electrons, 

A = 2\/2G F N e E v = 1.52 x 10" 4 eV 2 Y e p( g/cm 3 )£( GeV). (2.8) 

For antineutrinos, A is replaced with —A, and U with U*. Here Y e is the electron fraction and 
p(t) is the matter density. For neutrino trajectories through the earth's crust, the density 
is typically of order 3 g/cm 3 , and Y e ~ 0.5. For propagation through matter of constant 
density, the transition probabilities can be written in the form Eq. (|2.5|) if the mass-squared 
differences and mixing angles are replaced by the corresponding "matter" counterparts. Long 
baseline neutrino experiments are sensitive to matter effects, and the magnitude of the effect 
strongly depends on the baseline length and neutrino energy. Some examples are depicted 
in Figs. [I] and 0. 
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Figure 1: Transition probabilities for neutrinos (green, top curve) and antineutrinos (blue, 
bottom curve) in matter and vacuum (red, middle curve) as function of the energy for L= 735 
km, Am 2 3 = 3 • 1CT 3 eV 2 (normal hierarchy), atm = 7r/4, Am| = 1 x 10~ 4 eV 2 , #0 = 7r/6, 
|[L,| 2 = 0.04, and 5 = 0. 



There are some unknowns related to the neutrino mass pattern which can be addressed 
with the "help" of the matter effects. As alluded to before, the current data leave us with 
two alternatives for the spectrum of the three active neutrino species: a "normal" neutrino 
mass hierarchy or an "inverted" neutrino mass hierarchy. In the case of a "normal" mass 



hierarchy, the "solar pair" of states is lighter than u 3 , 



i.e. m 



2 > m 2 ,m 2 . In the case of 
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Figure 2: Transition probabilities for neutrino (green, top curve) and antineutrinos (blue, 
bottom curve) in matter and vacuum (red, middle curve) for 2 GeV as function of the 
distance, Am^ = 3 • eV 2 (normal hierarchy), # atm = 7r/4, Am^ = 1 x 10~ 4 eV 2 , 

e Q = vr/6, \U e3 \ 2 = 0.04, and 5 = 0. 

inverted hierarchy, the states of the solar pair are heavier than i/ 3 , i.e. m\ < m\ ~ m\. The 
key difference between these two hierarchies is then that, in the normal hierarchy, the small 
U e 3 admixture of v e is in the heaviest state whereas in the inverted hierarchy, this admixture 
is in the lightest state. The difference between both schemes is parameterized by the sign of 
Am^U A positive (negative) Am^ will point towards a normal (inverted) hierarchy. 

In going from v to V, there are matter-induced CP- and CPT- odd effects associated with 
the change A — > —A. The additional change U —>■ U* introduces further effects (this is the 
"genuine" CP- violation), which are usually sub- leading. Note that the matter effects depend 
on the interference between the different flavors and on the relative sign between A and Am 2 3 . 
Consequently, the experimental distinction between the propagation of v and V (the sign of 
A) can possibly determine the sign of Amf 3 . Fig. |3] shows the transition probabilities for the 
two possible choices of the sign of Am^ in the case of |Am 2 3 | = 3 • 10~ 3 eV 2 , a t m = tt/4, 
Am| = 1 x 10" 7 eV 2 , 9 Q = vr/6, \U e3 \ 2 = 0.01, and 6 = 0. 

What is currently known about the oscillation parameters? The atmospheric neutrino 

* Another way of treating the neutrino mass hierarchy is by defining ra\ > > m\, and redefining the 
solar, atmospheric and reactor angle depending on whether Amf 2 is larger or smaller than Am^. In such 
a scheme, the reactor data would limit |C43| 2 (normal hierarchy) or |f/ e i| 2 (inverted hierarchy). 
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Figure 3: — > v e oscillation probability for Am^ > (solid line) and Am^ < (dashed 
line). |Amf 3 | = 3 • 10" 3 eV 2 , # atm = vr/4, Am| = 1 x KT 7 eV 2 , 6 Q = vr/6, \U e3 \ 2 = 0.01, and 
5 = 0. 



puzzle points to a "well determined" mass-squared difference band (12 



1.5 • 10" 3 eV 2 < Am 2 tm < 5 ■ 10~ 3 eV 2 



(2.9) 



with sin 2 (26 l atm ) > 0.88, while the solar neutrino puzzle points to many disconnected regions. 
Loosely speaking, Am 2 & can lie anywhere between 10~ 10 eV 2 and 7 x 10~ 4 eV 2 (the best 
upper limit is provided by the CHOOZ experiment). Global fits to all data have historically 
identified the following well-known regions fl3H : 



LMA, large mixing angle: Am| = (1 - 10) • 10~ 5 eV 2 , tan 2 Q = 0.2-5; 



SMA, small mixing angle: Am| = (4 - 10) • 10 -6 eV 2 , tan 2 6 Q 



10) • 10 



-3. 



LOW-QVO, large mixing angle, small Am|: Am| = (0.5 — 2) • 10 7 eV 2 , tan 2 6 Q = 
0.1 - 10; 

VO, large mixing, vacuum oscillations: Am?, = (4 — 6)-10~ 10 eV 2 and (6 — 8)-10~ n eV 2 
tan 2 6 Q = 0.1-3. 



Notice that all the solutions satisfy 



Am% < Am 2 t 



(2.10) 



implying that there is indeed a hierarchy between the two distinct mass-squared differences. 
Current global fits point towards electron neutrino oscillations to active neutrinos (either 
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muon or tau neutrinos) in the LMA region (a small region in the LOW area is also allowed 
while the SMA and VO regions are strongly disfavored). Although the LMA solution is 
the most favored, this ambiguity must be removed by identifying the correct solution to the 
solar neutrino problem. This issue will be addressed by the KamLAND JE4[] reactor neutrino 
experiment, which is capable of confirming or excluding, unambiguously, the LMA solution, 
and precisely measuring the solar oscillation parameter if LMA is indeed correct. However, 
it is important to remember that if the mass-squared difference driving solar neutrino oscil- 
lations lies in the "high-end" of the LMA region, the length of the baseline will be too long 
and the oscillation peaks will not be resolved experimentally 115, 16|. The positron spectrum 



will be depleted (proportionally to the solar mixing angle) but will not be distorted as the 
oscillatory term will average out. In this case, the mass difference will be confined to a region 
whose lower limit is given by the KamLAND sensitivity and whose upper limit is dictated 
by the CHOOZ bound. Thus, as this mass-squared difference is an essential input in analyz- 
ing long baseline experiments, the capabilities of any future determination should take this 
possibility into account. Furthermore, a solar mass difference belonging to the upper part 
of the LMA region will have a strong impact on the oscillation probabilities, as the often 
neglected sub-dominant oscillations become not that sub-dominant any more. 

In order to determine the sensitivity of a given experiment to some parameter, it is 
imperative to distinguish different scenarios regarding the solar mass difference. We decide 
to study three different scenarios: 

1. KamLAND does not see an oscillation signal, implying that the LOW solution would 
be correct [15| | (or, perhaps one of the disfavored non LMA solutions). In this case, 



solar driven oscillations become negligible and there is no room for measuring CP 
violation in the neutrino sector. Nonetheless, one can still attempt to determine |t/ e 3| 2 , 
and, if a nonzero \U e ^\ 2 is observed, determine the neutrino mass hierarchy. 

2. KamLAND does provide a precise measurement of the solar mass-squared difference. 
The determination of the CP violating phase and \U e ^\ 2 from a simultaneous fit to both 
parameters becomes possible. 

3. KamLAND sees an overall suppression of the total rate but is not capable of measuring 
the mass-squared difference. In this case one must attempt to measure not only |£7 e 3| 2 
and S, but also the solar mass-squared difference. While the overall picture is rather 
"dirty," there is no reason to believe that a combined analysis is impossible (after all, 
for such large Ar/ig, there will be no shortage of z/ e -induced events!). 



§This is true only if CPT invariance holds. If CPT is broken KamLAND, using reactor antineutrinos, 
will not be able to constrain the solar (neutrino) spectrum (Tt| . 
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3 NuMI Off- Axis Beams 



The Neutrinos at the Main Injector (NuMI) [|TH tertiary beamline was designed to provide 



an intense beam to the MINOS experiment The z/^i^J's are derived from secondary 
7r + (7r~) and K beams that are allowed to decay within a 675 m decay tunnel. 120 GeV 
protons will be extracted from the Main Injector via a single turn extraction (8.6 /is pulse, 
cycle time 1.9 s.) and focused downward by 58 mrad, where the proton impinge into a 0.94 m 
graphite target to produce the secondary hadron (7r, K) beams. NuMI is expected to receive 
4 x 10 13 protons/pulse. 

MINOS is in a cavern of the Soudan mine located at a distance of 735 km from FNAL. 
The beam is tuned to make the experiment well aligned with respect to the beam axis. 
Another way of efficiently using this beamline would be to construct detectors located away 
from the beam axis. The resulting neutrino beam energy spectra at the different locations 
can be predicted from energy and momentum conservation in the ir decay process: 

E = rnl-ml = 0.004 GeV 

2(E W - cos 9 V ) - Ptt cos 9 V ) ' 

where m n and m M are the rest masses, E n and p n are the pion energy and momentum, and 
cos6 u is the angle at which the neutrino is emitted with respect to the pion direction. The 
maximum angle in the lab frame relative to the pion direction is related to the neutrino 
energy by: 

Or' = (30 + t T>MeV , (3-12) 

where 30 MeV is the neutrino momentum in the rest frame of the pion, while At takes 
into account the nonzero transverse momentum of the decaying it. As shown in Fig. f|(a), 
if 9 V ~ the neutrino energy is proportional to the pion energy (E v = 0A4E n ), while at an 
off-axis location (9 V ^ 0) there is a maximum neutrino energy which is independent of the 
energy of the parent pion. Therefore, the off-axis configuration allows one to use a fraction 
of the "total" beam that is characterized by having lower E v . The maximum flux for a fixed 
E v will be obtained when operating close to the corresponding 9™ ax , see Fig. |](b). The 
lower energy neutrinos provided by NuMI off-axis beams are highly desirable because they 
allow beams which are more suitable for oscillation studies, given the current knowledge 
of oscillation parameters (see Table [I]), while still having large enough energy and baseline 
length to be sensitive to matter effects (see Sec. ||). 

Two toroidal magnetic horns sign and momentum-select the secondary beam. The horns 
are movable allowing one to obtain different neutrino energy spectra. For example, Fig. |5] 
depicts the expected energy spectrum at different location for the low energy (horns 10 m 
apart) and the medium energy (horns 27 m apart) horn configurations. As shown, the off- 
axis beams are characterized by having a narrow and well-defined energy distribution with 
fluxes higher than the corresponding on-axis energy. In addition, the harmful high energy tail 
of the on-axis beams is not present, as expected from energy and momentum conservation, 
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Figure 4: top - allowed E v values for several 9 V values, bottom - shows 1/ cos#* that is 



proportional to the neutrino flux. cos#* = y'l — tan 2 0j, is the v production angle in the 
rest frame of the 7T. 
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Table 1: The first transition for — > vx occurs when 1.27|Amf 3 |(L/^) = n/2. Note 
that baselines much greater than 900 km are not allowed in our set up, because they would 
require "levitating" detectors. 



Am 2 3 (eV 2 ) 


0.002 


0.0025 


0.003 


Energy (GeV) 


1.5 (2) 


1.5 (2) 


1.5 (2) 


Length (km) 


928 (1236) 


742 (989) 


618 (824) 



Eq. ( [Up . All these distribut ions are calculated using; the GNuMI GEANT based Monte 
Carlo PDfl , and include the full beamline, target and decay pipe description. 

In addition to the off-axis beam properties described in Fig. [|, there are two basic facts 
that should be kept in mind when optimizing the beam configuration: (1) the beam flux 
decreases as the baseline increases, and (2) the beam flux decreases as the off-axis angle 
increases. Here we have performed two types of detector location optimization. 

The first one is referred to as the naive beam optimization, where we optimize the detector 
location such that 1.27|Am 2 3 |L/ E u = tt/2. That is, for a given |Am 2 3 | and L we have set 
8 exp ~ 9™ ax £ Q corresponding E u to be the maximum energy, which by definition 

will also correspond to the point of maximum flux. For example, 30 MeV/2.25 GeV ~ 
12 km/900 km ~ 0.013 mrad and 30 MeV/1.75 GeV ~ 12.5 km/735 km ~ 0.017 mrad. 
These results are depicted in Figs. |] and [7| for Am 2 3 = 0.003 and 0.0025 eV 2 , respectively. 

From this exercise we learn that, if we are interested in configurations with small Am 2 3 , 
very small energies will be required if the baseline is short. As a consequence, large angles 
are required, say < 20 mrad, and the lower acceptance in the medium energy beam for 
pions with wide angles becomes a limitation. This can be compensated for by increasing 
the baseline, therefore increasing the required energy and lowering the required angle. More 
important is the fact that, with the medium energy beam, there is a higher overall flux 
between 1.5 and 3 GeV, and a lower high energy tail than with a low energy beam, if one 
operates with angles that are smaller than ~15 mrad. These two facts make the medium 
energy beam preferable. 

The off-axis reduction of the high energy tail of the medium energy beam is depicted in 
Fig. [7|. There are two reasons for this: (1) the transverse momentum of the pions is smaller 
for the medium energy configuration making /\™ edlum < A 1 ™, and (2) the mean energy of 
the kaons in the medium energy beam is higher producing neutrinos that are of higher energy 
and therefore less harmful to the analysis, see Fig. [5]. 

In order to determine how well different transitions can be measured, all possible 
backgrounds have to be taken into account. Due to the short beam pulse duration and 
well defined repetition rate of the resulting neutrino beamline, it is possible to reject the 
cosmic ray background even if the detector considered here is located at the surface. On 
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Figure 5: On and off-axis beams for the low and medium energy NuMI horns configuration. 
A full beam simulation was made using the GEANT based GNuMI Monte Carlo. 
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Figure 6: Beam optimization for a fixed value of L/E v given Am 2 3 ~ 0.003 eV 2 . The 
medium energy configuration gives a higher event yield and cleaner beam conditions. 
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Figure 7: Beam optimization for a fixed value of L/E u given Am^ ~ 0.0025 eV 2 . The 
low (medium) energy configuration gives a higher event yield at 735(900) km. 
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the other hand, beam neutrinos induce other background events, which must be dealt with. 
The flavor composition of the "neutrino beam" is depicted in Fig. || For the neutral current 
(NC) background calculations, all neutrino flavors contribute. For — > ^-transitions, the 
inherent v e in the beam contributes with further (irreducible) background events. Assuming 
that we can not detect the sign of the outgoing lepton in charged current (CC) events, the 
relevant quantity is the ratio of [y e + i/ e )/(z/ M + z/ M ), which is depicted in Fig. [H]. The v T CC 



cross section for E u < 5 GeV is zero ||21|| , and therefore not relevant for the off-axis beam 
under consideration. 

Studies of matter effects and CP violation require the comparison of and oscillations. 
An "anti-neutrino beam" can be produced by reversing the polarity of the horns. This beam's 
flavor composition and [y e + u e )/{u^ + z/ M ) are also depicted in Figs. |9] and [U]. The beam 
related background is at the same level for both horn polarities, but the CC event rate is 
at least three times smaller for with energies above 2.5 GeV and five times smaller for 
the lower neutrino energies. As shown in Fig. [U], this is mostly caused by the cross section 
difference, and not by the tc + , iy~ production rate. 

In the second type of beam optimization, we look at locations with smaller 9 V . These 
locations are characterized by having a broader energy spectrum and a higher overall flux for 
neutrinos with energies below 3 GeV. This is illustrated in Fig. |B|(top), where we show the 
energy distributions at several radii and at a fixed length. Clearly, for a detector located at 
a radius of 10 km (13 mrad) instead of 12.5 km (17 mrad), we have at least two times more 
flux between 1.5 and 3 GeV. These broader energy beams are suitable for a wider range of 
Amf 3 , but how useful they are for — > v e appearance experiments will depend on how well 
one can keep the NC background under control. Therefore, in order to conclude how much 
better these larger 9 V off-axis experiments really are requires a full evaluation of the signal 
and backgrounds with a realistic detector simulation and reconstruction. We have performed 
such analysis and the details are given in Sec. |J Here we will just summarize the results by 
evaluating a figure-of-merit (FOM) at each location for — > v e in the case of full mixing 
and \U eZ \ 2 = 0.01. We defined the FOM as S/VS + B, where S and B are the signal and 



background events that survive all the cuts in the reconstruction. As depicted in Fig. [12 
off- axis experiments with angles between 10 < 9 V < 13 mrad from the axis have a high FOM 
for all values of Am^. The high FOM is not only due to the characteristics of the beam 
and oscillation probabilities (see Fig. |13[), but also due to the fact that in all cases we can 
keep the NC background at the 0.5% level, while the reconstruction efficiency is about 40%. 
If we look in detail at the case of Am 2 3 = 0.003 eV 2 , we can see that the naive beam tune 
performed at 735 km cannot compete with smaller 9 V locations at the same baseline. This is 
not true for the naive beam tune performed at 900 km, where at that location 9 V is already 
small enough to give us a high integrated flux. We can still obtain a 20% increase in the 
FOM by reducing the baseline to 735 km and 9 V to 10 mrad, but at the moment this is not 
particularly relevant, given the current lack of knowledge of neutrino oscillations parameters. 
For example, in Fig. [13] we have assumed a normal neutrino mass hierarchy, and an inverted 
hierarchy would significantly modify what the optimal conditions are. Furthermore, the 
inclusion of (potentially significant) solar oscillations and CP-violating effects further cloud 
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Figure 8: Contribution from ir and kaons to the beam composition for off-axis beams in 
the medium and low energy beam horn configuration. 
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Figure 9: Beam composition for positive (+) and negative (-) horn currents, left: for the 
a low energy configuration at L — 735 km and R = 10 km. right: for the a medium energy 
configuration at L = 900 km and R = 12 km. 



the picture. Nonetheless, the FOM is flat enough that any "reasonable" choice of baseline 
and opening angle should be "close" to optimal. We, therefore, will do all our analyses for 
a baseline of 900 km and at a radius of 11.5 km. 



3.1 NuMI Off- Axis Beams With A New Proton Driver 

The Proton Driver design described in [^, will allow us to bring the NuMI neutrino 
beam power up from 0.4 MW to 1.6 MW. This design is based on an 8 GeV circular machine 
with a circumference of 473.2 m, and it will provide 2 x 10 13 protons per pulse instead of the 
assumed 5 x 10 12 of the current booster. In addition, the total luminosity could be further 
increased by 30% if the current linac gets a 200 MeV upgrade. In this case, we would get 
3 x 10 13 protons per pulse. This machine is estimated to cost US$160M. 

An alternative design made out of only a linac to accelerate protons up to 8 GeV, using 



SNS and Tesla style superconducting cavities, is also under consideration [p5|. 



4 Detector Simulation and Expected Signal Efficien- 
cies 

To determine our — > v e detection capabilities, we are considering a highly segmented iron- 
scintillator detector. A full study of the expected detector performance, event reconstruction 
efficiencies, and background contamination was performed. 

The detector is made up of 4.5 mm thick iron foils (one quarter of one radiation length) 
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Figure 10: The {y e + v^jiy^ + v^) does not degrade when going from (+) to (-) polarity in 
the horns. In both cases, the v e + v e fraction of the beam is below 0.5% in the signal region, 
top - is for a low energy configuration at L = 735 km and R = 10 km. bottom: is for a 
medium energy configuration at L — 900 km and R = 12 km. 
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Figure 11: Comparison between the neutrino and anti-neutrino beams for charge current 
events. As shown, the lower rate is due to the cross section and not to the ir production 
rate. 



18 



0.9 



0.8 



0.7 



0.6 



0.5 



0.4 



0.3 



a L=735km, R=8km * L=900km, R=10km 

L=735km, R=10km 
▼ L=735km, R=l 1.5km 
o L=735km, R= 12.5km 



L=900km, R= 11.5km 
A 



A 



T 

O 



A 



A 



A 



O 



O 



▼ 

o 



o 



0.0025 



0.003 



0.002 



0.00275 



0.004 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 



Figure 12: FOM for — > v e appearance for different values of Am^ in the case of full 
mixing and |£7 e 3| 2 — 0.01 for different baselines and 9 V , but with the same detector and 
reconstruction criteria. We only select events with visible energy between 1 and 3 GeV. If 
this highly segmented detector were at the MINOS location and Am^ were 0.003 eV 2 , the 
corresponding FOM would be 0.54. This is to be compared to the FOM for the MINOS 
detector that is only 0.39 [22]. This means, that if we ignore the background uncertainty, 
about half of the gain in sensitivity is due to the location and the other half to the detector. 
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Figure 13: Transition probability for — > v e in the case of maximal atmospheric mixing 
and | t/ e 3 1 2 = 0.01, for different values of Am^. We assume Am^ 2 <C 10~ 5 eV 2 . 

interleaved with 1 cm thick scintillator planes. The scintillator strips are oriented at plus or 
minus 45° from vertical, alternating every other plane. The width of each readout cell is 2 cm. 
A 3 cm long air gap is left between each two iron-scintillator pairs to further improve the 
detector performance by increasing the effective radiation length at no additional cost. This 
choice is a compromise between the need for a good separation of two close electromagnetic 
showers on the one hand, and good clustering of each individual shower on the other. 

Signatures of v e and CC events, as well as of NC events, were studied in GEANT-based 
Monte Carlo simulations, with the GMINOS program developed at Fermilab [pq| . Typically, 
a 1-2 GeV v e induced shower will leave hits in 10-20 consecutive planes while ones will do it 
in around 40 consecutive planes, making possible a complete track finding procedure. High 
transverse segmentation provides good separation of two close tracks; preliminary studies 
showed that approximately 66% of 7r°'s at 1-2 GeV can be successfully identified as producing 
two well separated showers in at least one view. 

The event reconstruction consists of track fitting and track selection. Tracks are fitted 
and examined in each view separately. A good track is required to give hits in at least 4 
planes and have good \ 2 f° r a straight line. Most tracks coming from charged pions and 
recoiling protons can be rejected by requiring a mean track width of at least two cells and a 
width at maximum of at least three cells. "Baby tracks" , that is, tracks found in the vicinity 
of the end of a longer main track and less than half of its length, come mostly from secondary 
particles within the same shower and are discarded. Conversely, two tracks of comparable 
length and/or pointing at the same interaction vertex are a signature of a NC event with a 
7T° in the final state. 

After track selection, a signal candidate event is required to leave exactly one good track 
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Table 2: Expected event yields in the case of — > u e appearance for 1 kton-year for a 

detector located at different baselines (L) and distance from the beam axis (R), and different 

horn polarity (±). We assume that we are in vacuum, and maximal mixing solutions with 
A m 2 = Qm3 eV 2 and \ v _i2 = _ 01 _ 



L, R, ± 


735 km km + 




Signal CC v e 


0.709/3.351 = 0.212 




CC v e 


0.297/6.339 = 0.047 




CC u„ 


0.053/378.391 = 0.0001 




CC u T 


0.059/4.505 = 0.013 




NC 


0.613/165.686 = 0.0037 




L, R, ± 


735km 10km + 


735km 12.5km + 


Signal CC v e 


0.687/1.578 = 0.435 


0.342/0.897 = 0.381 


CC v e 


0.145/1.413 = 0.102 


0.105/0.825 = 0.127 


CC ^ 


0.012/26.97 = 0.0005 


0.0062/15.554= 0.0004 


NC 


0.132/32.1816 = 0.0041 


0.0411/18.363 = 0.0022 


L, R, ± 


900 km 11.5 km + 


900 km 11.5 km - 


Signal CC v e 


0.554/1.282 = 0.432 


0.183/0.460 = 0.399 


CC v e 


0.10145/1.036 = 0.098 


0.038/0.561 = 0.067 


CC 


0.0069/18.04 = 0.0004 


0.0027/10.924= 0.0002 


NC 


0.113/24.950 = 0.0046 


0.045/10.346 = 0.0044 



in each view. Additional selection criteria are imposed on the event basis. To optimize the 
ratio of v e signal to intrinsic u e background, a window in the total visible energy is defined 
(for Am^ tm = 0.0025 — 0.003 GeV 2 a reasonable choice is 1-3 GeV). A small missing pt 
with respect to the beam direction is required, a minimum fraction of the total event energy 
carried by the track (both criteria helping reject NC), and no track longer than 28 planes 
in any view (suggestive of a muon). The remaining NC background is further reduced by 
checking for a displacement of the beginning of the shower with respect to the interaction 
vertex, the latter being identified by the trace of the recoiling proton (if any) . 

The above criteria allow a rejection of about 99.7% of all NC events, while CC are 
suppressed to a negligible level. The resulting reconstruction efficiency and remaining con- 
taminations from NC and CC as a function of the incoming neutrino energy are shown 
in Fig. III. 



In order to evaluate how well we can do — > v e searches in the different off-axis beams 
locations, we have convoluted the beam spectra shown in Fig. ^| and the resulting recon- 
struction efficiencies shown in Fig. |14[ The results are summarized in Table ^| and Fig. [12 



21 



0> 0.2 
> 

° 

ox 

&0.002 
Cm 

© 

a 1 
© 

O.i 

to 0.05 








J 4 5 6 7 
v Charged Currents 



4 5 6 
Neutral Currents 



4 5 6 
v Charged Currents 



8 9 10 
E v (GeV) 




S 



iO 



E v (GeV) 



8 



10 



E v (GeV) 



Figure 14: Reconstruction efficiency for the selection of v e CC events, the NC background 
and the CC events for highly-segmented iron scintillator calorimeter. 
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5 Simulated Data Analysis 



5.1 Not the LMA Solution 

If KamLAND does not observe a suppression of the reactor antineutrino flux, the LMA 



solution to the solar neutrino puzzle will be excluded [27, 15 1, indicating that Am 2 2 <C 
10~ 5 eV 2 and/or tan 2 # <C 1. In this case, it is well known that the CP-odd phase 5 is 
not observable in standard long-baseline experiments, not only because solar oscillation do 
not have enough time to "turn on," but also because matter effects effectively prohibit any 
neutrino transition governed by the solar mass-squared difference. This being the case, one 
can only study v a <-> v$ transitions governed by the atmospheric mass-squared difference. 
For the reasons outlined earlier, we will concentrate on — > u e transitions and, possibly, 

-> V e - 

As mentioned in the Sec. 3, for the same running time, the number of interactions due 
to muon-type neutrinos obtained with the "neutrino beam" is significantly larger than the 
number of interactions due to muon-type antineutrinos obtained with the "antineutrino 
beam." Therefore, it seems logical to start running with the "neutrino beam" and decide 
whether one can observe an excess of z/ e -like events in the off-axis detector. This is done by 
simulating data and testing whether the observed number of events is significantly more than 
the the expected number of background events. We define the x _s ig ma average sensitivity 
region by 

, (DATA — BKG) 2 

where DATA is the averaged number of observed events for a given set of theoretical 
parameters, while BKG is the expected number of background events. cr BKG is the ex- 
pected error on the estimation of the number of background events. Unless otherwise noted 
a 2 BKG = (0.1 x BKG) 2 . In practice DATA = SIGNAL(U e3 , Am 2 , 5, ■ ■ -) + BKG. The "+1" 
accounts for the fact that we are computing the sensitivity of an average experiment (see 
|2§ for details). 

Fig. [T5| depicts the two and three sigma sensitivity to \U e -i\ 2 as a function of the number 
of kton-years of accumulated neutrino beam data collected off-axis, in the case Am 2 3 = 
±3 x 10~ 3 eV 2 , sin 2 6> at m = 1/2 and Am 2 2 = 10~ 7 eV 2 , sin 2 6 Q = 1/4 (for concreteness. The 
precise value of the solar parameters is irrelevant in the case at hand). The definition of 
one kton-year is the following: it is the amount of events observed after one year of running 
of the current NuMI beam (see Fig. |]) in a one kton detector at 900 km, 11.5 km off-axis. 
Many features are readily noticeable. First, in order to be sensitive to values of |£/ e 3| 2 which 
are significantly smaller than the current CHOOZ bound (|t/ e 3| 2 ^ 0.05 [^J), one is required 
to accumulate more than 40 kton-years of data, in the case of a normal hierarchy, or more 
than 150 kton-years in the case of an inverted hierarchy. This means, assuming the nominal 
NuMI beam, roughly two or eight years of running with a 20 kton detector. Note that the 
10% uncertainty on the background determination dictates that, even after accumulating 
an infinite amount of statistics, the three sigma reach of the off axis experiment plateaus 
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Figure 15: Two and three sigma sensitivity reach for \U e ^\ 2 as a function of the running 
time (in kton- years), for a normal neutrino hierarchy (solid lines) and an inverted neutrino 
hierarchy (dashed lines). |Am 2 3 | = 3 x 1CT 3 eV 2 , sin 2 6> atm = 1/2 and Am 2 2 = 1CT 7 eV 2 , 
sin 2 # = 1/4, 5 = 0. For a normal neutrino mass hierarchy, a 20 kton detector with 
(without) an proton driver upgrade will be able to detect a three sigma signal for z/ M — > v e 
oscillations, if \U eZ \ 2 < 0.0015(0.0028) after 5 years. 



at around \U e s\ 2 ~ 1 x 10~ 3 (2 x 10~ 3 ) for a normal (inverted) hierarchy. As is clear from 
Fig. in the case of an inverted hierarchy, the sensitivity is significantly worse. This is 
also expected, since matter effects enhance the u e appearance rate in the case of a normal 
hierarchy and reduce it in the case of an inverted hierarchy (see Fig. |3|). 

Note that the sensitivity would be significantly different for different values of |Am 2 3 | and 
that, by design, the sensitivity is optimal at around | Am 2 3 | ~ 3 x 10~ 3 eV 2 . We have verified, 
as discussed earlier, that it does not deteriorate significantly for |Am 2 3 | ~ (2 — 4) x 10~ 3 eV 2 . 

If one detects an excess of ^ e -like events, the next step is to determine the value of 
| t/ e 3 1 2 . Again, we do this by performing a \ 2 fit to the "data." We will assume that 
the atmospheric parameters |Am 2 3 | = 3 x 10~ 3 eV 2 , sin 2 # atm = 1/2 are precisely known. 
Fig. [L^(top, right) depicts x 2 — Xmin as a function of \U e s\ 2 corresponding to 120 kton-yearsQ 
of "data" collected with a neutrino beam (as defined earlier, the neutrino (antineutrino) 
beam consists predominantly of Note that, while the data were simulated with 

Am 2 3 = +3 x 10 -3 eV 2 and \U e -^ — 0.008, a different solution, with the same goodness of 
fit, is found for Am 2 3 = — 3x 10~ 3 eV 2 , |t/ e 3| 2 = 0.015.Q This implies that if the neutrino mass 

*This corresponds to six years of running with the current NuMI beam configuration and a 20 kton 
detector. With a proton driver, however, the same amount of data can be collected in 1.5 years. This will 
become crucial later. 

'It is important to reemphasize that |Ato^ 3 | is not a fit parameter. It is assumed to be known from 
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hierarchy is not known, instead of obtaining a precise measurement of \U e3 \ 2 = 0.008±0.0025 
(these are two sigma error bars), one is forced to quote a less precise (very non- Gaussian) 
measurement: 0.0055 < |t/ e3 | 2 < 0.018 at the two-sigma confidence level. 

The situation can be improved significantly if, after running with the neutrino beam, 
one runs with the antineutrino beam. Fig. |l6|(top, right) depicts x 2 ~ Xmin as a function of 
\U e s\ 2 corresponding to 300 kton- years of "data" collected with the antineutrino beam. As 
mentioned before, one is required to run much longer with the antineutrino beam in order to 
obtain a statistical significance comparable to the one obtained with the neutrino beam. One 
should readily note that 300 kton-years would correspond to 15 years (!) of running with the 
current NuMI beam configuration and a 20 kton off-axis detector. Here the presence of a 
proton driver becomes vital: it increases the beam intensity by a significant factor (nominally 
four), and reduces the 15 years to a bearable 3.74 years, such that the combined neutrino 
plus antineutrino beam time is slightly longer than 5 years. 

Again, the same behavior as before is observed: one obtains two distinct values of \U e3 \ 2 
depending on the assumption regarding the neutrino mass hierarchy, with a significant dif- 
ference: this time the effect is "reversed." The reason for this is simple: with the neutrino 
beam, the inverted hierarchy reduces the v e appearance signal compared to the normal hi- 
erarchy and, therefore, in order to correctly fit the data, a larger value of \U e z? (compared 
to the one obtained with the normal hierarchy) is preferred. In the case of the antineutrino 
beam, the inverted hierarchy enhances the v e appearance signal, and a smaller value of \Ue^ 
is preferred. This allows one to separate the two signs of Am 2 3 if the information obtained 
with both beams is combined. This is what is done in Fig. [l6|(bottom,left). Note that in 
this case the "wrong" model is about sixteen units of \ 2 away from the "right" model. It is 
also curious to note that, even with the wrong hypothesis, a similar measurement of \U e3 \ 2 
is obtained. This coincidence, which will not be considered too relevant, is a consequence of 
the fact that the data with the neutrino and antineutrino beams "pull" the measured \U e3 \ 2 
in opposite direction, and their combination meets somewhere "in the middle." 

Finally, in order to determine how well the two different signs of Amj 3 can be separated, 
Fig. |16|(bottom,right) depicts Xplus — Xminus as a function of the input value of |£/e 3 | 2 , plus 
the input Am 2 3 > 0. Note that for \U e3 \ 2 ^ 2x 10~ 3 , a y 2 separation of more than two units 
can be obtained. One can turn this into an exclusion of the "wrong" sign by noting that, for 
an average experiment, x 2 = 2 for the "correct" hierarchy, if one combines the data obtained 
with the two beams (2 is the number of degrees of freedom in this case). This implies that, 
for |£7 e3 | 2 = 0.01, the wrong hypothesis yields x 2 ~ 24, which is excluded at more than four 
sigma. A three sigma confidence level determination of the neutrino mass hierarchy would 
be obtained at \U e3 \ 2 ^ 0.005. 



In summary, if the LMA solution is excluded by KamLAND [27] (or, perhaps, a solar 
neutrino experiment, like Borexino [|3D|), future long baseline neutrino experiments can only 
hope to measure the magnitude of the \J e3 element of the neutrino mixing matrix, and deter- 
mine the neutrino mass hierarchy. On the other hand, both the estimated sensitivity reach 



different sources, possibly the study of the Un — * Vu disappearance channel in the off-axis experiment! 



25 




Figure 16: top-(x 2 — Xmm) as a function of | t^e3 1 2 , assuming both neutrino mass hierarchies, 
upon analyzing simulated data consistent with Am^ > and \U e ^\ 2 = 0.008 after 120 
kton-years of neutrino-beam running (left) and 300 kton-years of antineutrino-beam running 
(right), bottom, left-same as above, after combining the two data sets. bottom,left-difference 
of minimum value of x 2 obtained with the hypothesis Am^ > and Am^ < as a function 
of the input |t/ e3 | 2 . See text for details. |Am 2 3 | = 3 x 10~ 3 eV 2 , sin 2 atm = 1/2 and 
Am 2 2 = 10~ 7 eV 2 , sin 2 6 Q = 1/4, 5 = 0. 
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and the measurement of \U e 3\ 2 are very "clean," in the sense that they are not clouded by 
other physical effects (this will become clear after the next couple of sections). Furthermore, 
it is important to emphasize that both of these tasks are of the utmost importance, and are 
already enough to justify pursuing this kind of activity. 



5.2 LMA Solution with Spectral Distortions at KamLAND 



If the best fit point to the current solar data |T3| is close to reality, the KamLAND reactor 
neutrino experiment will be able to not only observe a depletion of the reactor antineutrino 
flux, but also determine the values of Am^ and sin 2 9 Q with very good precision \ Tf[ [3~5| , |31 



This being the case, in this section we address the issue of determining the neutrino 
mixing parameters \U e ^\ 2 and 8 at the off- axis experiment. As in the previous section, we 
start by determining the sensitivity of the off-axis experiment and 120 kton-years running of 
the neutrino beam to observing an excess of i/ e -like events. Fig. [IT] depicts the three sigma 
sensitivity in the (8 x |f/ e3 | 2 )-plane for Am 2 3 = 3 x 1CT 3 eV 2 | for 120 (300) kton-years of 
running with the (anti)neutrino beam. The sensitivity depends significantly on the CP-odd 
phase 8, and, as expected, the sensitivity is best for 8 ~ tt/2 in the case of running with a 
neutrino beam (8 ~ 37r/2 for the antineutrino beam), where the "interference" between the 
"CP-odd term" and the "|£/e3| 2 term" is constructive (i.e., one observes more events) and 
worse at 8 ~ 37r/2, where the "interference" is destructive. For smaller values of Am 2 2 , the 



'z-shape' and 's-shape' observed in Fig. [17j degenerate into vertical straight lines, such that 
the sensitivity will no longer depend on the CP-odd phase. 

If a signal is observed, one can attempt to determine the mixing parameters \U e s\ 2 and 
8. Similar to what was done in the previous section, the atmospheric parameters Am 2 3 = 
3 x 10~ 3 eV 2 , sin 2 # atm = 1/2 will be assumed known with infinite precision, and the same 
will now hold for the solar parameters Am 2 2 = 1 x 10~ 4 eV 2 , sin 2 9 & = 1/4. Furthermore, 
we will also assume that the neutrino mass hierarchy is known.0 This is done in order to 
not cloud the results presented here. Fig. |T8|(top,left) depicts the one, two, and three sigma 
measurement contours in the (|£7e3| 2 x <5)-plane obtained after 120 kton-years running with 
the neutrino beam. The simulated data are consistent with |?7 e 3| 2 = 0.017 and 8 = 7r/2. 
One can readily note that while \U e z\ 2 can be measured with reasonable precision, virtually 
nothing can be said about 8. Furthermore, the fact that 8 is not known implies that a 
measurement of \U e ^\ 2 irrespective of 8 is in fact less precise than what can be obtained if 
the solar parameters are not in the LMA region. 

In order to improve on this picture, it is imperative to prolong our "experiment" and take 



*As before, the sensitivity is close to optimal for Ato 2 3 = (2 — 4) x 10~ 3 eV 2 . One should keep in mind 
that if the atmospheric mass-squared difference turns out to be significantly different from this range, a 
different beam configuration has to be consider in order to optimize the appearance signal in the off-axis 
detector. 

§It may turn out, for example, that table top experiments [j§2| or the observation of supernova neutrinos 
p3| will be able to measure the neutrino mass hierarchy. 
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data with the antineutrino beam as well. Fig. |T8|(top, right) depicts the one, two, and three 
sigma measurement contours in the (\U e s\ 2 x <5)-plane obtained after 300 kton- years running 
with the antineutrino beam (as mentioned before, the longer running time is required in order 
to compensate for the "less efficient" antineutrino beam). Again, \U e3 \ 2 can be measured 
with some precision and nothing can be said about 5. A comparison of the two figures 
hints that a combined analyses may prove more fruitful. This is the case because while the 
neutrino beam yields a "z-shaped" measurement contour, the antineutrino beam yields an 
"s-shaped" contour. The reason for this is that the number of — > v e induced events is 
larger for 5 = tt/2 and smaller for 5 = 3tt/2. Therefore, the measurement will choose larger 
values of \U e3 \ 2 at around 5 = 3ir/2 in order to compensate for this small suppression. On 
the other hand, the number of —>■ u e induced events is smaller at 5 = tt/2 and larger at 
5 = 3n/2, and the opposite phenomenon is observed. 

Fig. |l8|(bottom,left) depicts the result of measuring \U e3 \ 2 and 5 using the combined 
neutrino and antineutrino beam "data." The situation is significantly improved, and now, 
a three sigma measurement of 5 ^ can be performed. Fig. |l8|(bottom,right) is similar to 
Fig. |I5|(bottom,left) , except that different input values of |f/ e3 | 2 ,(5 are chosen. As expected, 
the quality of the measurement is marginally worse for 5 = 3ir/2 (where S is consistent with 
zero at the three sigma level), and deteriorates as \U e3 \ 2 decreases. For 5 = 7r/2, one cannot 
determine that 5 ^ or it (i.e., no CP-violation) at the two sigma level if \U e3 \ 2 ^$ 0.004. It 
should always be kept in mind that the situation deteriorates for smaller values of Am 2 2 - 

It is worthwhile to note that an alternative (perhaps more standard) way to look at CP 
violation would be to construct an "asymmetry-like" parameter (e.g., something proportional 
to P(Vfj, — ► u e ) — P(i7^ — > v e )). This approach will not yield a better measurement of the 
oscillation parameters 8 and \U e3 \ 2 (as we are using all the experimental data and perform a 
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Figure 18: top - One, two, and three sigma measurement contours in the ( | t^e3 1 2 x ^)- 
plane, after 120 kton-years of neutrino-beam running (left) or 300 kton-years of antineutrino- 
beam running (right). The simulated data is consistent with \U e s\ 2 = 0.017 and 5 = tt/2. 
bottom, left - same as above, after the two data sets are combined. The solid star indicates 
the simulated input. bottom,right - same as before, for different simulated data points 
(indicated by the stars). Am 2 13 = +3 x 10" 3 eV 2 , sin 2 atm = 1/2, Am 2 12 = 1 x 10~ 4 eV 2 , 
sin 2 9 Q = 1/4. 
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"global fit"). Furthermore, it is not even an appropriate "direct detection" of CP- violation, 
given that our neutrino and antineutrino beams are not pure (especially the antineutrino 
beam) and matter effects contribute a significant amount of "fake" CP-violation. 

In summary, if KamLAND observes a suppressed and distorted reactor antineutrino 
spectrum, long baseline neutrino experiments can potentially measure both the magnitude 
of the U e 3 element of the neutrino mixing matrix and the CP-odd Dirac phase. Note that 
attempting to measure the CP-odd phase is not "optional" - the fact that it is unknown 
introduces a substantial uncertainty on determining It/^ 2 . In order to disentangle the two 
unknown parameters, it is crucial to take data using both a neutrino and an antineutrino 
beam.p] 



5.3 HLMA Solution, no Spectral Distortions at KamLAND 

If Am 2 2 ~ 2 x 10~ 4 eV 2 (note that this possibility is not currently excluded by solar neutrino 
data [|H|]), KamLAND will not be sensitive to the very rapid oscillatory pattern, and will 
only be able to observe an overall suppression of the solar neutrino flux [|15| , |T6|j . In this case, 
the mixing angle sin 2 9 & can be measured with some precision by determining the overall 
suppression factor, but the value of Am 2 2 will only be constrained to be larger than some 
lower limit. An upper limit will be provided by either future solar data or by the current 



CHOOZ data [29]. In order to be conservative, we will consider the latter case, and assume 
that AtUq <7x 10~ 4 eV 2 for large solar angles. 

If this scenario turns out to be correct, precise measurements of the atmospheric param- 
eters and the solar mixing angle will probably be available, while |t/ e 3| 2 , 5 and the precise 
value of Amj 2 will remain unknown.^] 

What are the consequences of having a very large but poorly measured Am 2 2 ? The 
biggest consequence, perhaps, is that even for very small values of |t/ e 3| 2 , a significant amount 
of zvlike events will be observed. This implies that the "sensitivity to |[/ e 3| 2 ," as discussed 
in the two previous sections is not a particularly meaningful quantity to study. Furthermore, 
as one may fear, this will also lead to a Am 2 2 versus \U e z\ 2 "confusion," (this was already 
alluded to in [Hi]) similar to the one observed between \U e 3\ 2 and 5 in the previous section 
(and which continues to exist here, of course). In other words, a moderate Am 2 2 and a 
large |£/ e 3| 2 will yield as many events as a large Am 2 2 and a small |C/ e 3| 2 . We already 
learned from the previous two sections that one will be required to run both the neutrino 
and the antineutrino beams (and accumulate enough statistics with both) in order to try 
to disentangle the three parameters. One should keep in mind that, while the situation is 
rather confusing, the number of observed events is going to be large for very large values of 
the solar mass-squared difference. 

^Another possibility would be to use two different detectors at different positions and/or a different beam 
with a different energy spectrum. We do not consider this possibility in this study. 

"A short-KamLAND or long-CHOOZ reactor experiment would certainly resolve this issue |$4|. Such an 
experiment has not been proposed yet (see, however, |35[] ). 
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In order to address what kind of measurement one may be able to perform under these 
conditions, we simulate data for 5 = tt/2, \U e3 \ 2 = 0.012, Am^ = 4 x 10~ 4 eV 2 , Am 2 3 = 
3 x 10~ 3 eV, sin 2 ^atm = 1/2, and sin 2 #0 = 1/4, assuming 120 kton- years of neutrino beam 
running and 300 kton-years of antineutrino beam running. As before, we assume during the 
data analysis that the atmospheric parameters, the neutrino mass hierarchy and the solar 
angle are known with infinite precision. 

The results of the three parameter fit are presented in Fig. [19], where we plot the three 
two-dimensional projections of the three sigma surface in the (Am 2 2 x |£7 e 3| 2 x<5)-space. Many 
comments are in order. First of all, one should note that the solar mass difference cannot 
be measured with any reasonable precision - it is only slightly better known than before, 
namely, it lies somewhere between 2 x 10~ 4 eV 2 (the KamLAND bound) and 6 x 10~ 4 eV 2 
(slightly better than the CHOOZ bound). The top left-hand panel depicts the Am 2 versus 
\U e3 \ 2 confusion alluded to earlier quite well. It is curious to note however, that the capability 
to decide whether \U e3 \ 2 is nonzero or not is not weak - the loose constraints on Am 2 2 are 
already enough to guarantee that \U e3 \ 2 > 4 x 10~ 3 . Most importantly, perhaps, at the 
three sigma level there is solid evidence that 5 7^ 0, n. The reason for this is that, for 
Am 2 2 = 4 x 10 -4 and 6 = tt/2, — > v e transitions are very suppressed, and a zero CP-odd 
phase would yield far too many i/ e -like events when the antineutrino beam is on. 

In summary, for very large values of Am 2 2 , KamLAND will not be able to measure the 
solar mass-squared difference with reasonable precision. Under these circumstances, long- 
baseline neutrino experiments are required to measure Am 2 2 along with |£7e3| 2 and 5. We 
find that, if \U e3 \ 2 is large enough, there is a reasonable chance that CP violation can be 
observed, while \U e3 \ 2 can be measured rather poorly. The measurement of Am 2 2 is even less 
precise. As before, accumulating enough statistics with both the neutrino and antineutrino 
beams is required, and, particularly on this case, one would profit immensely from other 
measurements with different distances and/or baselines (see, for example, [j3~6|). 



6 Summary and Conclusions 

The ultimate goal for the next generation of neutrino experiments will be to tie some of the 
"loose ends" of neutrino masses and mixings i.e., determine the neutrino mass hierarchy, 
measure (or further constraint) the "connecting mixing angle" (|£7 e 3| 2 ), and explore leptonic 
CP violation (if the solution to the solar neutrino puzzle lies in the LMA region). 

Candidate next generation experiments include shooting an intense conventional muon- 
type neutrino beam towards a detector which is located conveniently away from the main 
beam direction. "Off-axis beams" have several advantages with respect to their on-axis 
counterparts. They are more intense, narrower and lower energy beams, "void" of a large 
high energy tail, and provide clean access to z/ e -appearance, given the existence of a suitable 
detector with good electron identification capabilities. In order to fully take advantage of 
such tools, however, we have argued that two different beams (a predominantly and a 
predominantly beam) are essential. "Antineutrino beams," however, provide an extra 
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Figure 19: Projections of the three sigma measurement surface in the (|t/ e 3| 2 x 5 x Am^)- 
space, after 120 kton-years of neutrino-beam running and 300 kton-years of antineutrino- 
beam running. The simulated data is consistent with \U e s\ 2 = 0.012, 5 = tt/2 and Am^ = 
4 x 10~ 4 eV 2 . Am 2 3 = +3 x 10" 3 eV 2 , sin 2 # atm = 1/2, sin 2 ^ = 1/4. 
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experimental challenge: the antineutrino cross sections are suppressed with respect to the 
neutrino ones, such that, in order to obtain a statistically comparable antineutrino and 
neutrino data set, a sizeable amount of running time is required. The solution we find to 
this issue is a more intense proton source, along with a large (but realistically sized) detector. 

Specifically, we studied how well experiments in the NuMI beamline with a proton driver 
upgrade plus an off-axis detector can address these issues. In order to properly estimate 
the experimental response (e.g., signal efficiency as well as beam and detector-induced back- 
grounds), we performed a realistic simulation of the NuMI beam plus the response of a highly 
segmented iron detector, followed by a detailed "data" analysis. The combination of a new 
proton driver and the off-axis detector using the NuMI beamline can, in five years, improve 
considerably our knowledge of the neutrino sector. 

To properly assess the capabilities of such a set up, it is crucial to explore all the different 
physics scenarios, which will be (hopefully) distinguished by the current KamLAND reactor 
experiment. For different values of the solar mass-squared difference, we obtain different 
results for a five year program with an upgraded NuMI beam and a 900 km long baseline 
off-axis experiment: 

1 . KamLAND does not observe a suppression of the reactor neutrino flux - | U e s \ 2 can be 
measured with very good precision and the neutrino mass pattern can be established, 
as long as \U e s\ 2 few x 10~ 3 . We emphasize that even in this "less fortunate scenario" 
key information regarding neutrino physics will be obtained. 

2. KamLAND sees a distortion of the reactor neutrino spectrum - one should be capable 
of measuring \U e 3\ 2 with good precision and obtaining a rather strong hint for CP 
violation, as long as \U e s\ 2 ^ few x 10~ 3 , 5 close to either tt/2 or Zir/2. As a simplifying 
assumption, we assumed that the neutrino mass hierarchy would be already determined 
by other means. 

3. KamLAND sees an oscillation signal but is not able to measure Amg - one should 
be capable of measuring \U e ^\ 2 with some precision and obtaining a strong hint for 
CP violation as long as \U e ^\ 2 ^ 10~ 2 , 5 close to either tt/2 or 3ti/2, even if Awig is 
poorly known. AttLq cannot be measured with any reasonable precision. Again, as a 
simplifying assumption, we assumed that the neutrino mass hierarchy would be already 
determined by other means. 

It is also important to stress that, as discussed here in some detail, order 2 GeV neutrinos 
and a 900 km baseline are appropriate to measure fake as well as genuine CP violation (i.e., 
both the neutrino mass hierarchy and the Dirac phase 5 of the neutrino mixing matrix). 
Furthermore, we have argued that similar results should be obtained for different baselines 
(as long as the off-axis distance is appropriately modified) and for different values of the 
atmospheric mass-squared difference (see Fig. [12]). In particular, the latter implies that one 
need not wait until a very precise measurement of Am 2 tm is made in order to decide where 
the off-axis detector is to be located. 
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We conclude by comparing some of the results obtained here with similar studies per- 



formed for a future JHF to Kamioka neutrino beam 36 . A neutrino beam from the JHF 



proton source [36] aimed at the Super-Kamiokande detector [12] is a neutrino project with 



a 295 km baseline aimed to start in 2007-2008. Similar to the NuMI-off-axis project, the 
physics goals are to measure with an order of magnitude better precision (compared to 
MINOS and CNGS estimates |IJ) the atmospheric parameters (5(Am^ m ) < 10~ 4 eV 2 and 
5(sin 2 20 atm ) < 0.01), confirm <-> z/ r -oscillations or discover sterile neutrinos by mea- 
suring the neutral current event rate, and to improve by a factor of 20 the sensitivity to 
Vfj, — > z/ e -appearance. After a five year program, the JHF-Kamioka program should be 
able to exclude, at the 90% confidence level, — > v e transitions for |?7 e 3| 2 > 0.0015, 
while at NuMI with a 20 Kton off-axis detector we will exclude, at the two sigma confi- 
dence level, |t/ e 3| 2 > 0.00085 (0.0015) with (without) an upgrade proton driver (assuming 
AitLq <C 10~ 4 eV 2 and a normal neutrino mass hierarchy). The main difference between 
the two programs should come from the longer baseline proposed here, which allows the 
NuMI off-axis experiment (but not the JHF-Kamioka program) to cleanly try to address the 
neutrino mass hierarchy. 
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